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Different computational methods are employed to evaluate elastic (rotationally summed) 
integral and differential cross sections for low energy (below about 10 eV) positron scatter- 
ing off gas-phase C2H2 molecules. The computations are carried out at the static and static- 
plus-polarization levels for describing the interaction forces and the correlation-polarization 
contributions are found to be an essential component for the correct description of low- 
energy cross section behavior. The local model potentials derived from density functional 
theory (DFT) and from the distributed positron model (DPM) are found to produce very 
high-quality agreement with existing measurements. On the other hand, the less satisfactory 
agreement between the R-matrix (RM) results and measured data shows the effects of the 
slow convergence rate of configuration-interaction (CI) expansion methods with respect to 
the size of the Cl-expansion. To contrast the positron scattering findings, results for electron- 
C2H2 integral and differential cross sections, calculated with both a DFT model potential 
and the R-matrix method, are compared and analysed around the shape resonance energy 
region and found to produce better internal agreement. 
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I. INTRODUCTION 

The increase of interest in high-quality measurements involving antimatter has attracted,in re- 
cent years, the attention of experimentalists and theoreticians in the field of molecular physics 
[lli I2I bl y, cO. This has generated a wealth of new information on the nanoscopic behavior of 
a broad variety of molecular systems when they are made to interact with beams of positrons at 
thermal and near-thermal energies. To understand the interaction of positron beams with matter, it 
becomes important to also be able to distinguish to what extent the additional features of positron 
interaction with molecules (e.g. Ps formation and positron annihilation) are related to positron 
dynamics and to positron-electron correlation features. The study of even the simplest of such 
observables, e.g. the elastic scattering integral and differential cross sections occurring below the 
thresholds of Positroniun (Ps) formation in polyatomic gases, already provides a very useful test- 
ing ground for the theoretical and computational models which are currently employed to analyze 
positron-matter dynamics [2J. 

While the electrostatic interaction can in principle be described exactly by an essentially re- 
pulsive potential due to the molecular network of (electrons-i-nuclei), different approximations for 
correlation-polarization effects - the Vpcp potential - play an essential role in deciding the quality 
of the adopted theoretical model over the whole range of relevant distances between target and the 
positron projectile. As the projectile nears the target, in fact, the repulsive Coulombic core further 
slows it down while the attraction from the bound electrons increases and strongly modifies its 
motion in the intermediate range of distances via a correlation mechanism reminiscent of multiple 



scattering effects [^. This short-range effect should therefore be energy-dependent and nonlo- 
cal and would asymptotically give rise to charge-induced polarization effects, the leading term of 
which will be given via the dipole-polarisability of the target molecule 16]. The evaluation of the 
Vpcp contribution to positron-molecule interaction is therefore central to theoretical scattering cal- 
culations and its correct evaluation within cross section modelling studies is one of the stumbling 
blocks to the quantitative interpretation of existing experimental findings. 

In the present paper we therefore carry out a detailed comparison of the results of different 
theoretical approximations for the correlation-polarization forces by using the acetylene molecule 
as a benchmark systen, in view of the availability of good quality experimental data on this system 
and of its relatively simple structure as a polyatomic target. 

The structure of the paper shall be the following: in the next Section II we will report in some 
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detail an outline of the methods employed to generate correlation-polarization and static potentials 
while Section III will give and discuss our results for the elastic (rotationally summed) integral 
cross sections. The differential cross sections will be given by section IV and our conclusions will 
be summarized in Section V. 

II. SCATTERING EQUATIONS AND INTERACTION FORCES 

In order to carry out our comparison between different treatments of correlation-polarization 
forces, we have tested three different approaches: we have employed the Single-Center-Expansion 
(SCE) treatment of the scattering problem and included the Vpcp potential in two different ways, 
i.e. the Density functional modelling (DFT), used by us before [7^, and the Distributed Positron 
model (DPM) also introduced earlier on [8] and employed within the SCE treatment. We have 
then tested a Configuration Interaction (CI) procedure and implemented it within the R-matrix 
(RM) treatment of the scattering process [9]. Comparison with calculations by other groups jiol] 
will also be given for the sake of completeness. In the following we shall provide a short outline 
of each of the above computational methods. 

A. The SCE Scattering Equations 

In order to obtain the scattering cross sections for polyatomic molecules, we need to solve the 
Schrodinger equation for the total system 

{H-E)'^ = (1) 

at the total energy E, for the corresponding wavetunction ^. Here H is the total Hamiltonian given 
by 

H = Hmoi + K + V (2) 

where Hmoi, K and V represent the operators of the molecular Hamiltonian, kinetic energy of the 
scattered positron and interaction between the incident positron and the target molecule, respec- 
tively. The Hmoi further consists, in general, of the rotational and vibrational parts 

Hmol = Hrot + Hyih (3) 
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whereby we exclude, at the collision energies considered, both electronic excitations and the Ps 
formation channels. 

The total wavetunction \1' described in the molecular frame (MF) reference system in which the 
z axis is taken along the direction of the main molecular axis, is expanded around a single-centre 
(SCE) as 

^(ri...r2, rp|R) = ^^oj(ri...r^|R)y?(rp|R) (4) 

where 



(^(rplR) = Y: r;Vir,\R)X]:tir,) (5) 

lTTfJ,h 

In equation (4), ri represents the position vector of the ith electron among the Z bound electrons 
of the target, taken from the center of mass. The quantity '^moi represents the electronic wavefunc- 
tion for the molecular target at the nuclear geometry R. The continuum function (y9(rp|R) refers to 
the wavefunction of the scattered positron under the full action of the field created by the molec- 
ular electrons and by their response to the impinging positron as described below. Each uj/^ is 
the radial part of the wavetunction for the incident particle and the Xj^j^ are the symmetry-adapted 
angular basis functions discussed earlier \ which we will not repeat further here. 

The label vr stands for the irreducible representation (IR), /i distinguishes the components of 
the basis for each IR and {tt/}, respectively. 

Since the molecular rotations and vibrations are often slow when compared with the velocity 
of the impinging positrons considered in the present study, we may apply the fixed-nuclei (FN) 



approximation 



1211 that ignores the molecular term of Hmoi in equation (2) and fixes the values of 
all R at their equilibrium locations in each of the target molecules. We then solve the Schrodinger 
equation in the FN approximation, make use of the MF system rather than the space-frame (SF) 
reference system: the two systems are related through a frame transformation scheme given, for 
example, by [I2I . 

After substituting equation (4) into (1) under the FN approximation, we obtain a set of coupled 
differential equations for ui^ where, for simplicity, v represents (nfih) collectively: 

if-2 - + ^'1 ""M^) = 2Y.{lv\V\l'v')u,,ir,\R) (6) 

I "'^p ) I'v' 

with 
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(^^;|V|/'z;') = J drpXUrp) * y(rp|R)X,,,,(fp) (7) 

Solving equation (6) under the boundary conditions that the asymptotic form of ul?^, is repre- 
sented by a sum containing the incident plane wave of the projectile and the outgoing spherical 
wave we obtain the corresponding S-matrix elements, Sly^,. The actual numerical procedure we 



have employed to solve that equation was given in detail in [[13 



After transforming the MF quantities into the SF frame, the integral cross section (ICS) for the 
elastic scattering, rotationally summed over molecular rotations, is given by 



Q 



\ ^ \ ^ I rplv I 2 

2 Z^z^ l-'ivl 

Iv I'v' 



(8) 



B. The DFT Modelling of Correlation-Polarization 

The present treatment of the short-range part of the full Ypcp interaction was first applied to 
positron scattering problems by some of the authors (vl] and is based on constructing the correla- 
tion energy e^~P of a localized positron in an electron gas, and in further connecting it with the 
correct asymptotic form of the spherical dipole polarizability component contained within the full 
potential reported on the rhs of equation (6). The quantity e^~P had been originally derived by 



Arponen and Pajanne 11511 by assuming that the incoming positron can be treated as a charged im- 
purity at a fixed distance Vp in an homogeneous electron gas, which is in turn described as a set of 
interacting bosons representing the collective excitations within the random phase approximation. 
Based on their work, Boronski and Nieminen I16I1 have given the interpolation formulae of e'^^^' 
over the entire range of the density parameter r^, which satisfies the relationship ^iTr^p(r) = 1, 
with p(r) being the density of the Z bound electrons. 

The relationship between the correlation potential Ycorr and e^~'^, which is consistent with the 
local density approximation and a variational principle for a total collision system with the size of 
the target, is given by the functional derivative of e^^^ with respect to the electron density in the 
target: 

V,orr{rp\R) = ^{e''-'^[p{r,\R)]} (9) 
dp 
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where in our treatment /9(re|R) denotes the undistorted electronic density of the target, obtained 
from accurate Hartree-Fock (HF) calculations, as a a function of the molecular geometry and the 
re coordinates of the bound electrons. This quantity provides the probability of finding any of the 
molecular electrons near the impinging positron once an analytic ansatz is provided for e'^"^' |l61. 
Thus, the total Vpcp potential for the e+-molecule system can be assembled by writing, for each 
molecular geometry. 

V,,p{T,\R) = \ ^ ' ' " (10) 

[ Vpoi{rp\R) for Tp > 

The Vcorr is connected to the spherical part of the asymptotic form of the polarisation at the 
position of Tc, usually around 4.0 a^ for our systems. It corresponds to the outer crossing between 
the potential contributions of eq.(lO). 

The total interaction potential Vtot is therefore given by the exact static interaction Vgt between 
the impinging positron and the components (electrons and nuclei) of the molecular target , (for its 
detailed form see, for example, [|7D) plus the Vpcp given by equation [10] for each choice of fixed 



molecular geometry R: 

Vtot{rp\R) = Vst{rp\K) + Vpcp{rp\R) (11) 

C. The Distributed Positron Model 

An alternative model correlation-polarization potential we have used to describe the Vpcp in 



Eq. (10) is the distributed positron model (DPM) potential, K^f^^ [Il7l,ll8t]. The form adopted 



pep 

here for the correlation-polarization Vj^p^^ is based on a modification of the adiabatic polarization 
approximation, which provides a variational estimate of the polarization potential. In the adiabatic 
polarization approximation, in fact, the positron is treated as an additional "nucleus" (a point 
charge of -i-l) fixed at location Vp with respect to the center of mass of the atomic or molecular 
target. The target orbitals are allowed to relax in the presence of this fixed additional charge and 
the decrease in energy due to the distortion is recorded. The difference between the final and initial 
values of the energy defines the adiabatic polarization potential at one point in space. In order to 
describe carefully the spatial dependence of this interaction, the calculations need to be performed 
on a rather large number of three dimensional grid points. 
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However, due to nonadiabatic and short-range correlation effects, e.g., virtual Ps formation, 
the adiabatic approximation mayoverestimate the strength of the polarization potential for smaller 
values of r^, where the positron has penetrated the target electronic cloud. The present model 
corrects for this by treating the positron as a "smeared out" charge distribution rather than as a 
point charge. If the scattering particle really were an additional point charge, then the dominant 
short-range correlation effect would be virtual hydrogen atom formation into ground and excited 
states, and the 5-function distribution of positive charge at the center of mass would be correct. 
But, for a Ps atom, the positive charge is not localized at the center of mass, and to mimic this 
effect in computing the polarization potential we represent the positron as a spherically symmetric 
distribution of charge. This leads to a polarization potential that more closely reflects the correct 
physics and that smoothly reduces to the expected result for larger values of Vp without any need 
to select a crossover distance as done in the previous treatment. 

Within the above model, one can, in principle, choose any reasonable distribution that approx- 
imates the positive charge for virtual Ps embedded in the near-target environment. Studies by 
Gibson nivll have shown that constructing the positron charge distribution from the Is ST0-3G 



basis function with the tighter Slater exponent of ^=\.2A as recommended [|i8 



m for 



a molecular environment leads to accurate results. Once the V^^^^ potential is calculated, it is 
combined with the static potential to yield the total local interaction potential of Eq.(lO). 

After the development of the DPM to account for nonadiabatic polarization effects in positron- 



molecule scattering, a somewhat similar scheme was proposed by Bouferguene et al [|22l] for low 



energy electron-H2 collisions in which the polarization interaction is computed by replacing the 
impinging electron with a spherical Gaussian distribution of charge -1 . Like the DPM for positron 
scattering, this has the effect of reducing the overestimation of the adiabatic potential near the 
target and can be very efficiently implemented within a quantum chemistry framework. However, 
for electron scattering these authors ll22ll found it necessary to use a distribution that varies with the 
distance of the scattering electron from the molecular center of mass and involves a semiempirical 
parameter. 

In contrast, we have been able to obtain good agreement with various measured positron- 
molecule scattering results within the DPM using essentially the same procedure (i.e., without 



19, Q] as 



needing to adjust for details of the target molecule) on a variety of systems jlTl 
diverse as H2 and SFg. This past success is the chief motivation for including the DPM in the 
current study. In conclusion we construct the charge distribution from a Is ST0-3G basis function 
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with Slater exponent ^= 1.24 and once the V^^^ potential is calculated, it is combined with the 



static potential to yield the total local interaction potential of Eq. (10), as mentioned before ||23|]. 



D. The ab-initio R-matrix approach for e+/e scattering 



In the R-matrix method the space is divided into two regions: an inner region, defined by a 
sphere typically of radius 10 to 15 ao, and an outer region. In the inner region the complicated 
many-particle problem with correlation and polarization effects has to be solved. In the outer 
region the target is represented by a multipole expansion and the one-particle problem is solved 
by propagating the R-matrix outwards. The R-matrix provides the link between the two regions. 

For molecular targets, the calculations within the inner region reduce to a modified electronic 
structure calculation and standard quantum chemistry codes have been adapted to this purpose. In 
particular the UK molecular R-matrix codes use adapted versions [|28L 12911 of the Alchemy 
and Sweden-Molecule codes for diatomic and polyatomic targets respectively. All calculations 
discussed here employ updated versions of these codes. For diatomic targets an implementation for 



positron scattering was made some time ago [|3 
targets [jSJ. 



31(1 and has recently been extended to polyatomic 



The scattering wavefunction for a given energy E is built up as a linear combination CTII 



m{E) = Y.AK{E)mK (12) 

k 

where the coefficients are obtained by propagating the R-matrix into the outer region. The 
R-matrix basis functions are represented by the close-coupling expansion [|9|] 

= E E , (13) 

A i Br 

where the first sum runs over all products of target wavefunctions and electronic continuum 
orbitals ?7^, and A antisymmetrizes all electrons. The second sum runs over square-integrable 
functions $f of (Z + 1) particles. In the following we shortly describe the three different types 
of basis functions. The target wavefunction are obtained by diagonalizing the Hamiltonian 
for the target molecule containing electrons. In general the eigenfunctions of the target hamil- 
tonian are linear combinations of Slater determinants A^*^. In the scattering calculations the target 
functions are multiplied by the continuum orbitals r]i{rz+i) which are occupied by the scattered 
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particle 

A{e^%) = ^4i(A^'=(ri,---,r^)x77,(r^+i)) • (14) 

D 

Here the antisymmetrizer A acts only on the electrons. Therefore in electron-molecule scattering 
calculations the continuum orbital is anti-symmetrized with the target wavefunctions, whereas 
in positron-molecule scattering calculations the extra orbital is simply muliplied with the target 
wavefunction. In both cases the coefficients are kept fixed (see e.g. [33,1 for an efficient 
algorithm). The square-integrable functions are given by 

= <l>^f (ri,---,rz,rz+i) (15) 



where A^^ is a electron function, like a Slater determinant, and Xr(i"z+i) is a square-integrable 
spin-orbital. In the case of electron scattering the latter one is anti-symmetrized with the former 
orbitals. 

For electron scattering we have applied the static-exchange (SE) model and the static-exchange- 
plus-polarization (SEP) model. In both models the target is described by a Hartree-Fock wavefun- 
cion. In the SE model the additional electron can occupy all virtual orbitals as well as a set of single 
centered diffuse Gaussian orbitals which are used to represent the continuum within the R-matrix 
sphere 113 511 . In the SEP model all single excitations are added to the configurations generated in 
the SE model. 

In the case of positron scattering we use the same spatial orbitals for both electrons and 
positrons and define both a static and a static-plus-polarization (SP) model in analogy with the SE 
and SEP models employed for electron scattering. However, the positron can occupy all orbitals, 
including those orbitals, which are occupied by electrons. Furthermore, different spin-coupling 
rules apply in the two cases since the positron spin is not coupled with the bound electron spins. 



E. Semi-empirical R-matrix approach using an enhancement-factor 

In order to model correlation effects which are not fully described in the SP model, we have 
experimented with scaling the electron-positron attraction integrals by an empirically adjusted 
enhancement factor, /. These integrals are the ones which are routinely referred to as two-electron 
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integrals in standard quantum chemistry language. Here we wish to only increase the electron- 
positron attraction so do not alter the corresponding electron-electron integrals. 



{pq\rs)enh = f{pq\rs) (16) 
bp{ri)(l)g{ri) l--^]x^{rj)xl{r-)dridrj . (17) 



Here 0p(ri) and (j)q{ri) are electron orbitals, x^irj) and xII^t) are positron orbitals, and \r-^j\ = 
\vi — rj\ is the electron-positron distance. 

This form of scaling can be justified as follows. The second order contribution of a M0ller- 



Plesset type perturbative expansion of the electron-positron correlation energy Pq] is given by 



^(2) ^ ^ _ = E(^^l^)^'^ • (18) 



lata 



Here and ta denote energies of occupied and virtual electronic orbitals, respectively, ej and 
e-a are the same for positrons. Here we have introduced a second-order correction to the 

electron-positron attraction integral Since all denominators in the above expression are 

negative, this correction is always negative, and therefore has the same sign the integral 
itself. By defining the pair-dependent enhancement factor 



{ii 




[ii 


ii) 



the sum of first- and second-order contributions can be re-expressed as 

^(1+2) _ Y.i^i\n) + Y.{ii\nY^^ (20) 

a a 

a 

If we assume, that the enhancement is the same for each electron-positron pair, we can replace 

(2) 

the pair-dependent enhancement factor /.-r by an averaged enhancement factor / 



(21) 
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as used in this paper. In our computations we have used integrals of the more general type 

for which the second-order correction might be positive, resulting in enhancement factors smaller 

than one. We have however used the same enhancement factor for all integral types. 

For large distances between the positron and the scattering center the second-order contribu- 
tion to the electron - positron correlation energy goes over to the asymptotic polarization potential 



which is given by (in a.u.) BTIl 



for a linear molecule. Here ao (0^2) is the spherical (anisotropic) polarisability of the target 
molecule, and P2(cos 9) is a Legendre polynomial, where 9 is the angle between the vector linking 
the positron to the molecular centre-of-mass and that of the molecule. The long-range polarization 
is included automatically, if not completely, in calculations which use coupled-state expansions 



but not in the SP model. Below we also discuss the influence of including asymptotic 



polarization potential outside the R-matrix box, something that was also tested in earlier R-matrix 



studies of positron - molecule collisions pi . 



40(1 . In this paper the asymptotic polarization poten- 



tial has been neglected outside the R-matrix box. This effect has been studied in more detail in 



[|34ll . In all calculations the quadrupole moment generated by the target molecule was included in 
the outer-region. 

The present approach is related, but not equivalent, to scaling the positron charge or the elec- 
tron charge. Such an approach, in fact, would require a scaling of the positron - nuclear repulsion 
integrals or of the electron-nuclear attraction integrals in addition to the scaling done here. Fur- 
thermore, a scaling of the positron charge would require a scaling of the electrostatic interaction 
outside the box, whereas a scaling of the charge of the target electrons would introduce an at- 
tractive Coulomb-potential in the outer region. No such scaling of charges was done here and 
therefore no compensation of scaled charges is required. 

To use / = 1 products in the standard ab initio results of previous R-matrix studies of positron 
- molecule collisions. As shown below, the results are very sensitive to the choice of / and values 
of / only slightly bigger than unity yield surprisingly good agreement with existing experiments. 



13 



III. RESULTS AND DISCUSSION 

A. Low-energy positron scattering from C2H2 targets 

As mentioned in the Introduction, the present study intends to select gaseous C2H2 as a bench- 
mark system for an extended comparison of methods. Some of the earlier work on this molecule, 
both for positron and electron scattering studies, has been carried out by some of the present au- 



or POSH 



thors 11241 |25L |260. Therefore the present work employs the same molecular geometry to describe 
the equilibrium structure of the target, the same quality of Hartree-Fock basis set and, whenever 
necessary, the same value of the dipole polarisabilitv coefficients for the long-range polarization 



potential employed in the previous studies 



ies 



25 



The numerical convergence of the SCE has been carefully checked both on the multipolar 
expansion of the potential (lmaa;=30) and on the number of partial waves describing the scattering 
e+/e^ particle. 

In the R-matrix calculations for electron and positron scattering we have used the DZP basis 
set of Dunning and coworkers, an R-matrix box of a radius of 10 ao, and the continuum basis set 



of Faure et al., optimized for a box size of 12 oq [1350. We have taken the equilibrium geometry 
(rcc = 1.208A and TcH = 1.058 A) optimized with the the Hartree-Fock method using this basis 
set. For the positron-scattering calculations using the enhancement factor, we have empirically 
found a factor / = 1.004 by comparing the calculating integral cross section with the experimental 
results. Differential cross sections were calculated using the program POLYDCS H. 

We report in figure 1 the elastic (rotationally summed) integral cross sections for positron scat- 



ts with experiments from Sueoka and Mori 14211 The 



tering, comparing various theoretical resu 

earlier computations from Carvalho et al. 11^ are also given in the figure. The following consid- 
erations could be made from a perusal of the results reported: 

1 . The calculations at the static level, carried out with the same basis set and target description, 
provide the same results when using the SCE and the R-matrix treatment: both curves are 
given in the lower part of the figure and are essentially coincident, which is a good test for 
the two different scattering codes; 

2. When correlation-polarization forces are included, we see a dramatic change in the size 
and behavior of the cross sections obtained within the SCE approach and a remarkable 
agreement with the experiments down to very low collision energies; 
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3. The use of the DFT and DPN approximations for Vpcp potential yield results which are in 
very good agreement with each other and in good accord with the experiments. The values 
of the DPM cross sections are consistently smaller than the DFT data especially around 1 
eV and below. However, this difference keeps well below the experimental cross section in 
the considered energy range; 

4. The ab initio inclusion of polarisation effects within the R-matrix calculations using the SP 
model is seen, on the other hand, to be still far from convergence although it does give an 
upturn in the cross section dependance on collision energy near threshold. It is thus clear 
that the use of configuration interaction methods for treating correlation effects in positron 
scattering requires a much more extensive inclusion of additional states in order to be able 
to reproduce the experimental findings; 

5. The introduction of the empirical enhancement factor into the R-matrix method markedly 
improves the agreement with experiments of the calculated cross sections at all energies. 
Compared with DFT and DPM calculations the RM cross sections are now slightly larger 
at low energies and slightly smaller at higher energies, although they display the correct 
energy dependence. 

If we now turn to the angular distributions from elastic scattering processes, the results of figure 
2 report the behavior from the calculations using the DPM model within the SCE treatment of the 
scattering problem. The two panels of figure report the differential cross sections (DCS) over a 
broad range of collision energies; the lower panel shows their values at higher energies from 4.0 to 
10.0 eV, while the upper panel reports the same quantities down to 0.5 eV. The general trend of all 
angular distributions remains fairly similar at all energies and indicates the presence of an intensity 
"dip" that moves to lower angles as the energy increases. Thus, it is only a shallow minimum in 
the DCS values beyond 130° at 0.5 eV while it moves down to 30° and becomes more marked by 
the time the collision energy goes up to 10 eV. 

The DCS calculated using the R-matrix method via the SP-model, with and without enhance- 
ment factor, are given figure 3. The results using the SP-model without the enhancement factor 
are given by thick lines and do not show a strong forward peak. The calculations with an enhance- 
ment fact, as in the DPM calculations, show a strong forward peak. Indeed it would appear that at 
scattering energies above 1 eV the smaller cross sections given by the RM calculations are caused 
by their underestimation of the strong forward peak. 
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All computed cross sections turn out to be characterized by a large forward peak and appear to 
converge to a similar limiting value for all energies considered. In other words, the dominance of 
the 1=2 Legendre polynomial, which is associated with the dipole polarizability coefficient, shows 
up in the angular shape of the DCS at small angles and in the energy dependence of a sort of 
"magic" angle in the angular distributions. 

For the R-matrix calculations including the enhancement factor, the value / = 1.004, has been 
used as a result of optimizing the integral cross sections discussed before. The corresponding 
results in figure 3 are given by thin lines: as mentioned before, the inclusion of the enhancement 
factor introduces a much more marked forward peaks at all energies. Furthermore these DCS are 
also showing a minimum which moves to lower angles at increasing energy, as found by the DPM 
results of figure 2. At energies above 4 eV, however, a second minimum appears, which is not 
observed in the DPM calculations. 



Differential cross section for C2H2 have also been calculated by Carvalho et al [|43ll . and mea 



sured by Kauppila et al [|44|] . The angular distributions in both works are showing strong forward 
peaks, in agreement with both the DPM calculations of fig. 2 and the R-matrix calculations with 
the enhancement factor in figure 3. 

B. Electron scattering from gaseous C2H2 targets 

As mentioned before, one of the objectives of the present work was to revisit the quantum 
observables associated with scattering electrons for processes involving the same target, the C2H2 
molecule, in order to compare the performances of different as well as to investigate the differences 
with the previous comparison carried out for positron-C2H2 scattering attributes. 

The results reported by figure 4 therefore present the calculations at the exact Static+Exchange 
level carried out within the SCE expansion but treating both the Static and the exchange potentials 
as exact contributions since the latter was introduced as a discrete, converged expansion over a 



separable, additional basis set: for details see our discussion given in reference [|45ll . The dotted 
and dash-dotted curves reported by figure 4 clearly show the very good numerical agreement 
exhibited by the SCE and RM methods at the SE level of modelling: thus, these data further 
confirm the internal reliability of the two approaches in describing e+/e^ scattering processes. 

The inclusion of correlation-polarization forces, both within the SCE-DFT and RM approaches, 
is reported by the dashed and solid curves, respectively. Here again, the two treatments exhibit 
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reassuring agreement and similarities: the location of the open-channel (shape) resonance feature 
falls between 2.2 and 2.5 eV, the overall size of the cross sections is also rather similar over the 
whole range of energies and both curves fall out asymptotically to a very similar high-energy 
value. The threshold behaviour, on the other hand, exhibits marked differences very close to zero 
energies, although both methods would require further numerical tests within that region before 
anything could be decided on the origin of the discrepancies. 

Another interesting set of comparisons is reported by the data shown in the three panels of 
figure 5: we see there the behaviour of computed and experimental differential cross sections as 
a function of collision energies and for three different values of the scattering angles. The top 
panel refers to ^9=40°, the middle panel to ^=60° and the bottom one to 'd=90°. The solid lines 
refer to the R-matrix calculations while the dashed curves report the SCE-DFT calculations. The 
experimental data of Kochem et al. are from ref. [4Q. 

We clearly see again that the treatment of correlation effects for electron scattering processes 
are easier to model ab initio within the R-matrix approach than their treatment for positron scatter- 
ing attributes: the RM results are fairly similar to the DFT calculations without having to introduce 
any enhancement factor. Furthermore, both computational models follow reasonably well the ex- 
perimental findings, especially for the case of the DFT results while the RM results turn out to 
be slightly larger in the energy region of the shape resonance. Both methods, however, reproduce 
well the general energy dependence of the measured cross sections. 

IV. PRESENT CONCLUSIONS 

In the present work we have tried to analyze in some detail the relative performances of different 
computational treatments for obtaining low-energy scattering observables (integral cross sections 
and angular distributions in the elastic channels) associated with both positron and electron low- 
energy collisions with a specific polyatomic target: the C2H2 molecule at its equilibrium geometry. 

In particular, we have first made sure that both methods are comparable when the interaction 
forces are artificially simplified to be given either by Static interaction only (for e+ scattering) of 
by Static-i-Exchange interactions (for e~ scattering): the present results clearly show that indeed 
both sets of codes produce the same integral elastic (rotationally summed) cross sections in spite 
of their implementing very different computational procedures. 

The next step has been to add within both treatments the effect of correlation-polarisation forces 



17 



as implemented within each code. 

The results which compare our findings with experiments in the case of the positron scattering 
show that the DFT/DPM modelling of correlation forces, neither of which are entirely ab initio 
are both able to yield remarkable agreement with observed quantities, while the ab initio RM 
approach moves in the correct direction but manages to reach the same good agreement only when 
an empirical enhancement factor is included. Thus the RM approach still requires further work for 
positron scattering studies on the treatment of correlation-polarisation effects in order to achieve 
an acceptable level of convergence. 

The corresponding analysis of computed angular distributions largely confirms the above find- 
ings. The two methods produce very similar results only when the empirical enhancement factor 
is employed within the RM calculations. 

Finally, the comparison between scattering attributes, as obtained by the two methods analysed 
in the present work, has been extended to electron collisions from the same molecular target. The 
results reported by figures 4 and 5 indeed confirm the methods' reliability at the SE level but also 
indicate the quality of the RM outcome when electron projectiles are considered: both methods 
now yield similar results without using an enhancement factor within the RM calculations. In 
general, however, the DFT angular distributions appear to follow more closely the experimental 
data over a broad range of energies. 

In conclusion, the present benchmark calculations have allowed us to attain the following, 
interesting results: 

1. that the existing implementations of the R-matrix and of the SCE-DFT approaches for treat- 
ing e''"/e~ electronically elastic scattering off polyatomic targets provide essentially the same 
results when using the same interaction potentials; 

2. that the two codes model very differently the further inclusion of correlation-polarization 
forces and therefore show up at this level their differences of behaviour; 

3. that the DFT approach is found to be able to model such effects in ways which turn out to 
reproduce very well experimental results for positron and for electron scattering; 

4. that the RM approach shows faster CI convergence in the case of electron scattering data 
while still requires much larger expansion treatments in the case of positron scattering, as 
indicated by its empirical need of including an enhancement factor in the latter case. 
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We have therefore shown that at least two different approaches to multichannel scattering meth- 
ods which provide observables for polyatomic targets can be reliably employed in the future to 
yield complementary information for the ever increasing range of polyatomic gases that are being 
experimentally analyzed. 
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FIG. 1: Integral cross sections for positron-C2H2 scattering using various theoretical methods. The experi- 
mental results of Sueoka and Mori 14211 . and the calculations of Carvalho et al jlOll are also given. 
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FIG. 2: Computed e+ — C2H2 elastic differential cross sections using the DPM modelling of correlation 
forces, see text for details. 





FIG. 3: e"^ — C2H2 elastic differential cross sections calculated with the R-matrix method using the SP- 
model with (thin lines) and without (thick lines) using the enhancement factor (/enh = 1-004). 




FIG. 4: Electron scattering: computed symmetry components of the elastic integral cross sections (ICS) for 
electron scattering (upper panel) and a comparison of total ICS obtained either using the DFT correlation 
polarization model (soUd Une) or the static + exchange method only (dashes) 
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FIG. 5: Computed — C2H2 angular distributions over the location of the shape resonance energy and 
for three different angular values {9 = 40, 60, and 90 degrees from top to bottom panel). Experiments are 
from ref . 14611 



